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Abstract 

We analyze the influence of the characteristics of the STM-tip (applied voltage, tip radius) on 
the electroluminescence spectra from an STM-tip-induced quantum dot taking into account the 
many-body effects. We find that positions of electroluminescence peaks, attributed to the electron- 
hole recombination in the quantum dot, are very sensitive to the shape and size of the confinement 
potential as determined by the tip radius and the applied voltage. A critical value of the tip 
radius is found, at which the luminescence peak positions as a function of the tip radius manifest 
a transition from decreasing behavior for smaller radii to increasing behavior for larger radii. We 
find that this critical value of the tip radius is related to the confinement in the lateral and normal 
direction. 

PACS numbers: 73.21.La, 78.60.Fi, 73.21.Ac, 78.20.Bh 
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I. INTRODUCTION 



Scanning-tunnelling microscopy techniques have developed into an important tool for 
studying semiconductor nanostructures in recent years. The scanning-tunnelling microscope 
(STM) covers a wide field of applications. An example of such an application is the so- 
called tip-induced quantum dot, a special and unique type of quantum dot that can only 
be studied using a STM . When a negative bias is applied between the metallic 

STM-tip and the semiconductor sample, the electric field extends into the structure, and a 
hole accumulation layer is formed at the surface under the tip. Such an accumulation layer 
effectively screens the electric field. For sufficiently small tip apex radii, quantization occurs 
both in the radial direction and along the axis normal to the surface of the structure. It 
should be noted that the tip-induced band bending confines only a single type of carriers, 
holes in our case. The electrons are repelled from the region under the tip. By putting 
a barrier just below the surface, injected electrons are confined to the surface layer and 
radiative recombination with holes in the dot becomes possible. 

In this paper we analyze the influence of the characteristics of the confinement 
potential, determined by the STM-tip radius and the applied voltage, on the elec- 
troluminescence spectra of a tip-induced quantum dot in the GaAs layer of the 
GaAs/Al .25Ga .75As/In .25Gao.75As quantum well structure ^]. 



II. MODEL 



In our model we consider a multilayer structure sketched in Fig. 1. The shape of the 
STM tip is modelled by a truncated cone because, when the tip made of a soft metal (like 
Pt in our case) is pressed to the semiconductor layer, it becomes flat. The STM-tip and 
the GaAs layer (cap layer) are separated from each other by an insulator layer of width 
l\. The holes and electrons are confined in the cap layer of width l^. The width of the 
Alo.25Gao.75As barrier is The Alo.25Gao.75As layer borders on the Ino.25Gao.75As quantum 
well. The contact area of the STM-tip with the structure is a circle of radius i?tip- Owing 
to the axial symmetry of the system, cylindrical coordinates r, ip, z are used. The z-axis is 
directed along the symmetry axis of the tip, which is chosen parallel to the normal to the 
surface of the structure. While the hole motion is confined in all directions by the potential 
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well, the electron motion is confined only along the z-axis. 

In order to find energy spectra of holes, trapped in the quantum dot, we solve self 
consistently the Poisson equation 

Id ( dV (r, z)\ d 2 V(r,z) p h (r,z) 



r- 



(1) 



r dr \ dr J dz 2 e £ j 

where ph (r, z) = en h (r, z) is the hole charge density and £j is the dielectric constant of the 
material, together with the Schrodinger equation governing the hole motion: 

h 2 ( 1 1 d f d\ 1 Id 2 d 1 d\ T . , . 



2 y m ^J ( 2 ) r $ r V ^ r / m)j (z) r2 df 2 dz mf (z) dz 
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Here the index j labels the hole band type: j = 1 for a hole, which is heavy for the motion 
along the z-axis and light in the plane of the GaAs layer and j = 2 for a hole, which is 
light for the motion along the z-axis and heavy in the plane of the GaAs layer. For the case 
of a quantum well, the light and heavy holes exactly decouple at ky = 0. Since the lateral 
size of the quantum dot is much larger than the GaAs unit cell, only states around ky = 
contribute to the wavefunctions in the quantum dot. The effect of coupling between the 
valence bands is very weak in the close vicinity of the center of the Brillouin zone ky = 
(see, e.g., Refs. 0, |(|). The index s labels subbands due to the size quantization of the 
hole degree of freedom along the 2-axis. The index n is the radial quantum number and m 
is the angular quantum number of the lateral (in the plane of the GaAs layer) motion of 
the hole. The kinetic energy of the hole motion along the z-axis is larger than the kinetic 
energy of the lateral motion. (Our calculations show that the charge of holes is concentrated 
in a relatively thin (~ 3 nm) layer near the interface insulator/ GaAs. The radial size of 
the smallest dot considered in this paper is 12 nm.) Hence the adiabatic approach can 
be used, assuming that the hole degree of freedom along the z-axis is "fast" while the 
lateral one is "slow". Therefore, the hole wave function is represented in a product form: 
^j,s,n,m{r,4>,z) = ^f s (z-,r)^j snm (r)e im<l> , where ^f s (z;r) and ^\ s , n , rn {r)e im<t ' are the wave 
functions describing the motion along the z-axis and the lateral motion, respectively. The 
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Schrodinger equation governing \P+ a (z;r) reads 

where V(#, r) = et/(z, O+KxL-^) with U(z, r), the electrostatic potential. V^ r (z) represents 
band offsets for a hole |7[. 

We solve the Schrodinger equation © using the method described in Ref. 0, 0|. It 
yields the spectrum of the subband energies Ej iS (r) and wavefunctions ^fj- (z',r). Each of 
the energies Ej a (r) determines the top of a hole subband in the structure and plays the role 
of an adiabatic potential for the lateral motion. 

Consequently, the equation, which describes the lateral motion, is as follows: 

~2~Vt ( r ^) ~ ^) *^ (r) + Ej ' s (r) v "'-" ! " (r) = E ^^Ln, m (r)- ( 4 ) 

Aftarllg this nation a8 de 8 cnbe d to Re, Q, we obtain the ene rgy spectrum of the 
holes, confined in the STM-tip-induced quantum dot. 

Given a particular number of holes in the quantum dot, we find the quasi- Fermi level E F 
for holes from the equation: 

= £ E , E . 2 . E ^ • ( 5 ) 

3=1,2 s,n,m exp I ^ I + 1 

Then, within the Hartree approximation scheme, the hole density may be extracted from 

2 

n h (r,z)=J2 J2\yj,s,n,m(r,(l),z)\ 2 • (6) 

j=t* t^m exp F ) + 1 

This is the basic formula employed to invoke a self-consistent solution to the Schrodinger 
and Poisson equations. 

In order to include in our scheme the exchange effects we use the local density approxi- 



mation (LDA) as described in 
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where e^° m is the exchange energy of the homogeneous gas per particle. 

The exchange energy of a hole gas was evaluated in Ref. j 1 31 ] - The result for exchange 
energy per hole is 

E wx <■- 3 
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where the numerical function £ (w) is given by 



i 

rf + ! ^ + !to + lH(1 -^ Ha | + i(1 -*)/* Ha | 

cw = ^ (9) 

A parameter w represents the ratio between the light- and heavy-hole Fermi wave vectors: 



w= — = J — . 10 

&Fh V m h 

In the local density approximation, the exchange potential is related to the exchange energy 
e£T (n) by Q 

V«(r) = iA[ n ^»( n )]( r ). (11) 
The exchange potential is thus 

^x C (r^)=-— — (-^ ^|^- s , n , m (r,0,z)| 2 ) CM, (12) 

This formula is employed to include exchange corrections within LDA in the self-consistent 
solution of the Schrodinger and Poisson equations. When V(r, z) is substituted by V(r, z) + 
V xc (r,z), we need a straightforward numerical solver as shown in Fig. 2. The iteration 
procedure, involved to calculate the hole states and the corresponding electrostatic potential 
is subdivided into four steps. In the first step, the electrostatic potential V(r, z) in the 
whole structure is calculated in the absence of holes. Then the obtained potential is used 
as the initial guess in the iteration procedure. In the second step, the Schrodinger equation, 
describing the hole motion is solved. The obtained wave functions are used to calculate the 
hole density. In the third step, solving the Poisson equation with the charge density obtained 
in the previous step yields a new approximation for the electrostatic potential, which is 
involved in the next iteration. In the fourth step, we use the charge density obtained as a 
result of the second step to calculate the exchange potential Eq. ([12)1 . which is also used in 
the next iteration. At a given number of holes in the tip-induced quantum dot, the second, 
the third and the fourth steps are repeated consecutively. This sequence of steps continues till 
the maximal absolute value of the difference between the values of the electrostatic potential 
at consecutive iterations becomes less than a certain threshold (< 0.1 mV), which establishes 
the measure of the accuracy. When increasing the number of holes in the quantum dot, the 
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electrostatic potential, obtained earlier for a smaller number of holes, is used as an initial 
guess for calculations. This procedure guarantees a continuous link between the solutions 
obtained at the consecutive numbers of holes. 

Given the electrostatic potential formed by the hole charge, we solve the Schrodinger 
equation for electrons, 

h 2 f 1 1 d ( d \ 1 i d 2 d i d \ 

2 [m)l (z) rdr\ r dr) + m \\ {z) dtp 2 + dz (z) dz ) ' Z) 

+ V e (z,r)^ e Se ^ me (r, ( p,z) = E e Se ^ m q e Se ,n e ,mSr,(j),z), (13) 

where V e (z,r) = —eU(z,r) + V£ ar (z\ U(z,r) is the electrostatic potential, and V^ ar (z) 
describes band offsets for an electron (7|, 

6 U S m V [M ^ G ^ ^,if.GAl . 25 Ga . 75 As. 1 ] 

The index s e labels subbands due to the size quantization of the electron motion in the direc- 
tion normal to the surface of the structure. The index n e is the radial quantum number and 
m e is the angular quantum number of the lateral motion of the electron. The electron wave 
function is represented in a product form: ^ ne m (r, (p, z) = tyj-f(z; r)^f t n etm< ,{r)e ime ^ ^ 
where fyjrf(z; r) and ^sf,n e ,m e (r)e me ^ are the wave functions describing the motion along 
the z-axis and the lateral motion, respectively. 

Our numerical approach to solve the Schrodinger equation, which governs the lateral 
degree of freedom of the electron, is described in Ref. The set of energies E e s n m 

and wave functions \^ neme (r, <fi, z) together with the set of hole energies Ej tS>ntm and wave 
functions ^ j, s ,n,m(^, <t>, z ) are used to describe the electron- hole radiative recombination. The 
intensity of electroluminescence at a frequency Q in the quantum dot created by the STM 
tip is described by the following expression [5J]: 

I (HQ) ~ 22 fh(j,s,n,m) f e (s e ,n e ,m e ) P (j,s,n,m,s e ,n e ,m e ,hQ) , (15) 

j,s,n,m s e ,n e ,m e 

where fh(j,s,n,m) and f e (s e ,n e ,m e ) are, respectively, probabilities of finding holes and 
electrons in the corresponding quantum states, P (j, s, n, m, s e , n e , m e , HQ) is the probability 
of the electron-hole recombination. 
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We assume that the distribution functions fh (j, s, n, m) is described by the Fermi-Dirac 
function with the quasi-Fermi level E F : 



Following Ref. |9[, we assume that due to the radial accelerating field the non-equilibrium 
electron distribution is a strongly decreasing function of the angular quantum number m e . 
Hence, the electron distribution function is approximated as 



III. RESULTS AND DISCUSSION 

The characteristics of the structure under consideration are as follows. The thickness of 
the vacuum barrier is l\ = 0.5 nm, while the thickness of the GaAs and Alo.25Gao.75As layers 
are, correspondingly, l 2 = 17 nm and Z3 = 46 nm. The respective dielectric constants are 
e 3 = 1.0, e 2 = 13.2 and e 3 = 12.2. 

Figure 3 illustrates the position of the calculated peaks in the electroluminescence spec- 
trum as a function of the voltage V^ ip applied between the STM-tip and the semiconductor 
structure at 4.2 K. Clearly, the peak positions as a function of the bias between the sample 
and the tip reveal a red shift. This is explained as follows. An increase of the voltage shifts 
the bottom (top) of the dot (antidot) upwards nearly proportional to the applied voltage. 
Since the number of holes increases with the voltage, the width of the dot in the z-direction 
is reduced due to a stronger screening. Hence, an increase of V ti p leads to a deepening of 
the potential well for holes and to a rise of the barrier for electrons. This, in its turn, raises 
the energy of the electron states, which contribute to the recombination, with respect to 
the bottom of the conduction band in the absence of the applied voltage. The holes expe- 
rience both the rise of the energy level due to the shift of the bottom of the dot and the 
lowering of the energy level due to a stronger confinement. The net effect of both trends is 
an upward shift of the hole energy levels. The upward shift of the hole energy levels with 
increasing Vh d is larger than that for electrons. Thus, these shifts result in a red shift in the 

Q 

electroluminescence spectrum (see also Ref. 
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s e ,n e ,m e ) ~ S me0 . 
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In Fig. 4 we present the peak positions for the spectrum shown in the inset of Fig. 3 as 
a function of the radius of the STM-tip. When the radius of the STM-tip R t i P > 12 nm, one 
notices that the enlargement of the radius Z? t ip leads to a rise of the transition energies (blue 
shift). This is understood on the basis of the fact that an enlargement of the contact area 
shifts on average the conduction and valence bands upwards, resulting in a rising potential 
barrier for electrons and a slight deepening of the potential well for holes. Consequently, 
the recombining electron has to occupy a higher energy level in order to have a considerable 
overlap integral with a hole wave function. Our calculations show that the upward shift of 
the hole levels is smaller than the increase of the energy of the recombining electron. So, 
a change of the contact area influences both hole and electron levels, but due to the weak 
confinement for holes in the radial direction, the shifts of the electron and the hole levels 
jointly lead to a Uue *fi of the tuition ene rgy (this behavior of the tuition energy as 
a function of the tip radius has been analyzed in Ref. |9|]). 

As distinct from that type of behavior, for R tip < 12 nm an increase of the contact area 
between the STM-tip and the semiconductor sample results in a red shift of the transition 
energy. The explanation of this phenomenon is as follows. When the radius of the STM-tip 
i?ti P is much smaller than the thickness of the GaAs layer, the electrostatic field from the 
STM-tip hardly extends into this layer, as shown schematically in Fig. 5(a). Hence, the band 
bending in the GaAs layer is negligible. In this case the energies of the transitions between 
hole and electron levels are practically equal to those in the case without any external field. 
An enlargement of the contact area increases the extension of the electrostatic field into 
the structure, as shown in Fig. 5(b). This effect is similar to that of an increase of the 
applied voltage. In both cases transition energies become smaller (red shift). However, 
under the further enlargement of the contact area the electrostatic field within this area 
extends almost homogeneously over the whole semiconductor layer, see Fig. 5(c). As the 
luminescence originates in this case from the edge of the STM-tip- induced quantum dot, 
there occurs a smaller rate of decrease of the transition energies with increasing radius. 
Finally, when the radius i? t i P > 12 nm, the effects giving rise to the blue shift, described in 
the previous paragraph, prevail. 

Thus, we have shown a strong effect of the confinement potential, as determined by the 
STM-tip radius and the applied voltage, on the electroluminescence peak positions. Namely, 
while an increase of the voltage on the STM-tip results in a red shift of the electrolumines- 
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cence peaks |9( , an increase of the contact area between the STM-tip and the semiconductor 
structure results in a non-monotonous variation of peak positions. Below a certain critical 
value of the tip radius, the transition energies decrease with increasing radius, while above 
the critical radius of the tip, the transition energies are an increasing function of radius. 
These results demonstrate vast possibilities of tuning optical characteristics of the STM-tip 
induced quantum dot by varying parameters of the STM-tip. 
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FIGURE CAPTIONS 



Fig. 1: Scheme of the tip/insulator/GaAs/Al .25Gao.75As/In .25Gao.75As quantum well struc- 
ture. The shown hole density distribution is calculated for the STM-tip voltage —3.1 V 
and the radius of the tip 20 nm. 

Fig. 2: Flow diagram of the numerical calculations. 

Fig. 3: Positions of peaks in the electroluminescence spectra, calculated for the STM-tip- 
induced quantum dot at 4.2 K, R tip = 12 nm, versus bias V tip . Inset: Typical electro- 
luminescence spectrum calculated for the voltage on the STM-tip Kip = —3.1 V and 
the STM-tip radius R t i p = 12 nm. 

Fig. 4: Positions of the electroluminescence peaks from Fig. 3 as a function of the STM-tip 
radius at the voltage on the STM-tip V t i P = —3.1 V. 

Fig. 5: Sketch of the electrostatic potential profile for the STM-tip induced quantum dots 
with small (a), intermediate (b) and large (c) radius. Darker regions correspond to a 
larger absolute value of the electrostatic potential. 
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FIG. 1: M. D. Croitoru, V. N. Gladilin, V. M. Fomin, J. T. Devreese, M. Kemerink, P. M. Koenraad, 
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FIG. 4: M. D. Croitoru, V. N. Gladilin, V. M. Fomin, J. T. Devreese, M. Kemerink, P. M. Koenraad, 
J. H. Wolter 




FIG. 5: M. D. Croitoru, V. N. Gladilin, V. M. Fomin, J. T. Devreese, M. Kemerink, P. M. Koenraad, 
J. H. Wolter 
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